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We analyze the necrosis growth due to thermal coagulation induced by laser light absorption 
and limited by heat diffusion into the surrounding live tissue. The tissue is assumed to contain a 
certain tumor in the undamaged tissue whereof the blood perfusion rate does not change during 
the action. By contrast, the normal tissue responds strongly to increase in the tissue temperature 
and the blood perfusion rate can grow by tenfold. We study in detail the necrosis formation under 
conditions typical for a real course of thermal therapy treatment, the duration of the action is taken 
about 5 minutes when a necrosis domain of size about or above 1 cm is formed. In particular, if 
the tumor size is sufficiently large, it exceeds 1 cm, and the tissue response is not too delayed, the 
delay time does not exceed 1 min, then there are conditions under which the relative volume of the 
damaged normal tissue is small in comparison with the tumor volume after the tumor is coagulated 
totally. 



LASER INDUCED HEAT DIFFUSION 
LIMITED TISSUE COAGULATION 



Thermal coagulation of living tissue caused by local 
heating due to laser light absorption is one of the novel 
thermotherapy techniques of tumor treatment which is 
currently undergoing clinical trials (see, e.g., |l|]). Ther- 
mal coagulation is used to form a necrosis domain of 
desired form for the removal of the malignant tissue. 
So mathematical modeling of the necrosis growth is re- 
quired, first, to find out the physical limitations and the 
basic characteristics of the treatment and, second, to op- 
timize the therapy course. However, living tissue is ex- 
tremely complex in structure, thereby, for the adequate 
theoretical model to be developed and for the mathemat- 
ical modeling of the given process to be implemented reli- 
ably typical limit cases should be singled out and studied 
individually. 

One of these cases is the laser induced heat diffusion 
limited tissue coagulation [Q, || . It is characterized by the 
following features (Fig. Q). Absorption of laser light de- 
livered into a small internal region of living tissue causes 
the temperature to attain such high values (about or 
above 70 °C) that lead practically to immediate thermal 
coagulation in this region. Heat diffusion into the sur- 
rounding live tissue causes its further coagulation, giv- 
ing rise to the growth of the necrosis domain. In this 
case heat diffusion plays a significant role in the necro- 
sis growth because the necrosis size 5ft exceeds the depth 
of laser light penetration into the tissue. Therefore, the 
temperature distribution inevitably has to be substan- 
tially nonuniform and for the tissue to coagulate at pe- 
ripheral points heat diffusion should cause the tempera- 
ture to grow at these points. The latter property singles 
out the specific mode of thermal coagulation under dis- 
cussion from other possible types of thermotherapy treat- 
ment and that is why we refer to the necrosis growth 
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FIG. 1: The necrosis growth due to local thermal tissue coagu- 
lation limited by heat diffusion in the surrounding undamaged 
tissue, d = 1, 2 corresponds to the surface plane applicator 
and the internal cylindrical one. 



under the given conditions as to thermal tissue coagula- 
tion limited by heat diffusion. In particular, it turns out 
that the optimal implementation of this thermotherapy 
mode is characterized by the formation of necrosis do- 
mains of size 5ft > 1 cm and by the treatment duration of 
^course ~ / D ~ 2-8 min (where D w 2 • l(T 3 cm 2 /sec 
is the tissue temperature diffusivity) [§, || . 

Heat diffusion in the live tissue is affected substantially 
by blood perfusion causing the heat sink Q|. Thus, the 
characteristics of the spatial distribution and the dynam- 
ics of the blood perfusion rate should have a substantial 
effect on the necrosis growth limited by heat diffusion. 
Therefore, in modeling this therapy mode one has to 
take into account the tissue response to the temperature 
growth which can locally give rise to a tenfold increase in 
the blood perfusion rate B. The latter effect, in partic- 
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ular, is responsible for a substantially nonuniform distri- 
bution of the blood perfusion rate and visually manifests 
itself in a red ring ( "hyperemia ring" ) appearing around 
the necrosis region during the thermotherapy treatment. 

In the previous papers || [|, ^, |l0|, ll], [l^, |l3| basing 
on the free boundary description we have developed a 
model for the heat diffusion limited thermal tissue coag- 
ulation and studied the properties of the corresponding 
necrosis growth. In particular, we have shown that the 
given mode of thermal coagulation comprises two stages, 
fast and slow. At the former the necrosis domain grows 
fast and its size attains values of order 5ft~l cm. The 
latter is characterized by a substantially slower growth 
of the necrosis domain. Exactly at this stage the necro- 
sis growth is governed by the competition between the 
heat diffusion into the surrounding undamaged tissue and 
the heat dissipation caused by blood perfusion. There- 
fore, a sufficiently strong tissue response to temperature 
variations which, in turn, gives rise to an essential lo- 
cal increase in the blood perfusion rate has to affect the 
necrosis formation. Namely, the higher is the local blood 
perfusion rate, the smaller is the size of the necrosis do- 
main. 

Inside tumors the blood vessels are abnormal, in par- 
ticular, they have lost the ability to expand substantially 
in the response to tissue heating. As a result, inside a tu- 
mor the blood perfusion rate can increase no more then 
twice under strong local heating, although in the unaf- 
fected tissue the blood perfusion rate inside a tumor typ- 
ically exceeds the perfusion rate in the normal tissue || . 
So, if the normal tissue exhibits the strong response to 
temperature variations then the boundary of the necro- 
sis domain will penetrate deeper insider the tumor than 
in the normal tissue during the thermal therapy course. 
Therefore in this case, first, the relative volume of the 
damaged normal tissue can be less than the volume of 
the tumor after its coagulation. Second, when the necro- 
sis boundary reaches the tumor boundary the necrosis 
growth should slow down. At first glance it would seem 
that the necrosis growth is mainly confined to the tumor 
space. Exactly this phenomenon is the subject of the 
present paper and will be called the self-localization of 
the necrosis growth in active living tissue with a tumor. 



calized in a small region adjacent to the applicator whose 
thickness is substantially less than the characteristic size 
5ft > 1 cm of the necrosis domain formed during a typical 
course of thermotherapy treatment. 

In modelling heat propagation in the tissue we single 
out three regions: the necrosis domain Q n , where the 
blood perfusion rate is equal to zero 



j(r, t) = for r G 



(i) 



the undamaged tumor tissue Qt, and the undamaged 
(live) normal tissue Q\. Inside the necrosis domain Q n 
the tissue temperature T(r, i) obeys the heat diffusion 
equation for solids: 



9T - 2rp 



q, r e Q n . 
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where c and p are the specific heat and density of the tis- 
sue, and n is the cellular tissue heat conductivity. In the 
undamaged tissue (including also the undamaged part of 
the tumor) the temperature is governed by the bioheat 
equation (see also 

dT 



impera 
, 0): 
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Here T a is the temperature of arterial blood in systemic 
circulation, j v is the value of the blood perfusion rate j 
locally averaged over spatial scales 
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where L n = ln(//a) ^3 — 4 and l/a is the mean ratio 
of the individual length to radius of blood vessels form- 
ing peripheral circulation, the cofactor / ~ l A/ X^ ~ 0.5 
takes into account the counter-current effect ]l5|, |l6| , and 
the factor F — F{L n ) ~ 2 allows for the renormalization 
of the heat conductivity due to blood flow. It should 
be pointed out that the scale £ v (r) is exactly the mini- 
mal scale on which equation (||) is justified Q and the 
characteristic size of the necrosis domain formed dur- 
ing this mode of thermal therapy can be estimated as 

5ft > /QV 



II. FREE BOUNDARY MODEL OF THE 
NECROSIS GROWTH 
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(5) 



We study the necrosis formation caused by laser light 
absorption in living tissue with a tumor applying to the 
following model |, | | g 0, |ll|, g| 



The laser 

light absorption causes heat generation and, as a result, 
thermal coagulation of the tissue. For simplicity sake the 
heat generation rate q(r) is considered to be a beforehand 
given function of the spatial coordinate r that is indepen- 
dent of the tissue state and the time t (certainly, for t > 
and t = is the beginning of the action on the tissue). 
Besides, we assume the heat generation rate q(r) to be lo- 



Finding the relationship between the averaged and true 
blood perfusion rates, j v and j, we have to take into 
account that the scale £ v of averaging in its turn depends 
on the local value of j v (expression (Q)). This dependence 
enables us to specify this relationship in the form [|l4[ : 



jv ~ — V 2 lnj v 
cp 



3, 



r G 



Ql\jQt 



(6) 



where X v ~ is also a constant of order unity. 

Equation (0) should be subjected to a certain boundary 
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condition at the interface T of the necrosis domain be- 
cause it makes no sense to average the blood perfusion 
rate over the necrosis domain impermeable to blood. The 
physical layer separating the necrosis domain and the un- 
damaged tissue where the local vascular network is not 
damaged is complex in structure and contains a spatial 
increase of the blood perfusion rate from zero to the value 
in the undamaged tissue. In order to avoid the problem 
of analyzing the blood perfusion rate in this layer we take 
into account the following simplifying circumstance. On 
one hand, the typical size of the necrosis domain formed 
during a thermal therapy course and the characteristic 
length of temperature variations are of the same order 
about 1 cm. So, particular details of spatial variations in 
the blood perfusion rate on scale much less than 1 cm are 
not the factor. On the other hand, the damaged part of 
the vascular network located inside the necrosis domain is 
most probable to be made up of an artery and vein having 
supplied previously this region with blood as a whole and 
of shorter vessels formed by their branching. Therefore, 
the region containing the vascular network part in which 
blood flow is remarkably disturbed because of the necro- 
sis formation does not exceed substantially the necrosis 
domain. The latter enables us not to make difference be- 
tween the given layer and the interface T and to deal with 
a sharp jump of the blood perfusion rate at the necrosis 
interface. The desired boundary condition imposed on 
the averaged blood perfusion rate j v meets the require- 
ment that the normal gradient of the averaged blood per- 
fusion rate be equal to zero at the interface T: 



(7) 



We note that the adopted boundary condition will not 
hold if a large vessel goes through the necrosis domain. 
However, the probability of this event is small enough 
and this case should be analyzed individually. 

Now let us specify the tissue response to temperature 
variations. Only the undamaged normal tissue is consid- 
ered to respond to local heating, inside the undamaged 
tissue of the tumor the blood perfusion rate maintains 
the initial value jo assumed for simplicity sake to be the 
same in the normal tissue and the tumor 

j(r,t)=j for reQ t . (8) 

Blood vessels can expand only to a certain extent as the 
temperature grows. So when it becomes high enough, 
T > T vr , the blood perfusion rate j(r,i) attains a large 
but finite value j max and remains approximately con- 
stant. Keeping in mind this feature we describe the nor- 
mal tissue response to local heating by the equation: 

r^|+j$(T)=jo for reQi. (9) 

Here r is the delay time of the tissue response and the 
function $(T) is of the form 

] k=k for T<T - do) 

for T > T vr 



$(T) 



I a+(l 



where a = jo/jmax and the temperature at which the 
blood vessels exhaust their ability to expand is estimated 
as T vr « 45-46°C. 

Equations ®) and (||) governing the evolution of the 
tissue temperature should be completed by the boundary 
conditions relation the temperature field and the necrosis 
interface T. To solve this problem we have developed the 
free boundary model that assumes the temperature and 
the heat flux to have no jumps, i.e. the temperature 
distribution meets the following boundary conditions 
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and the normal velocity of the interface V is given by the 
expression 
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Here the constant ~ 1 and the function 

~T-T 



uj(T) = uoo exp 



(12) 



(13) 



is actually a convenient approximation of the Arrhe- 
nius dependence of the thermal coagulation rate oo(T) oc 
exp{ — ^}. The available experimental data [|l7j for the 
temperature dependence of the exposure time enable us 
to estimate the value of A as A - 3-5 °C (A ~ 3.6 °C for 
pig liver at To = 60 °C). Below in numerical calculations 
the dependence (|l3|) will be taken in the form 



;(T) = 0.2 exp 



T-60 



3.6 



mm 



(14) 



where the temperature T is in degrees Celsius. 



III. RESULTS OF NUMERICAL SIMULATION. 

SELF-LOCALIZATION OF THE NECROSIS 
GROWTH IN LIVING TISSUE WITH A TUMOR 

The stated model has been used to analyze numeri- 
cally the necrosis growth in the case shown in Fig. [|. 
The laser light absorption near the applicator gives rise 
to the heat generation rate q(y) uniform along the x-axis 
and decreasing with y as q(y) oc exp(— y/l&d), where the 
characteristic thickness Z a d of the laser light absorption 
layer is chosen less then the characteristic size 5ft > 1 cm 
of the necrosis domain formed during the thermal ther- 
apy mode under consideration, Z a d < 5ft. The intensity 
of laser light has taken such values that the tissue tem- 
perature directly near the applicator, y = 0, be about 
80 to 90 °C ( q(y) = 2exp(— y) for large tumor and 
q(y) = 4exp(— Ay) for small tumor). We have studied 
two opposite cases specified by the relation of the scale 
5ft (see expression (|J)) and the tumor size C. For a small 
tumor, £ < 5ft, the self-localization phenomenon cannot 
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FIG. 2: The system where the necrosis formation has been 
studied numerically based on the free boundary model. The 
xy-region of the plane and cylindrical symmetry imitates the 
effect of the surface applicator and that of the internal cylin- 
drical one. 



ficiently small tumor also cannot substantially affect the 
necrosis growth, as is seen in Fig. || exhibiting the necro- 
sis growth in living tissue with strong response without 
delay. The tumor radius was chosen equal to C = 0.5 cm 
the tissue response to heating was assumed to be strong, 
a = 0.1, and without delay, r = 0. This case has been 
considered because under such conditions the effect of 
large tumor on the necrosis growth is most pronounced. 
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come into being and this case is presented for illustra- 
tion only. For a large tumor, C > 5ft, there are condi- 
tions under which the self-localization is pronounced as 
is demonstrated directly in the present section. Besides, 
we have considered the ^-region of plane geometry as 
well as cylindrical one in order to compare the necrosis 
growth for applicators of different geometry, the surface 
plane applicator and internal cylindrical one. In the lat- 
ter case the radius ro of the applicator also was assumed 
to be small in comparison with the scale 5ft. 

In numerical simulation we have used the following 
typical values of the thermal conductivity k ~ 7 • 10 -3 
W/cm-K, the heat capacity c ~ 3.5 J/g- K, and the den- 
sity p ~ lg/cm 3 of the tissue, as well as set the blood 
perfusion rate jo ~ 0.3 min -1 and the factors / = 0.5, 
F = 2. Figure [| demonstrates the necrosis growth in 
living tissue without response. In this case the existence 
of the tumor has no effect on the necrosis growth. A suf- 



FIG. 3: The necrosis growth in living tissue without response, 
where the existence of the tumor has no effect on the necrosis 
formation. The dotted lines show the position of the necrosis 
interface at different moments of time. 



When the tumor size is sufficiently large, C > 5ft, and 
the tissue response is not too delayed, r <C 5ft 2 /D the 
self-localization phenomenon comes into being as demon- 
strated in Figures |-@. Namely, Fig. | exhibits the necro- 
sis growth in living tissue containing a tumor of size of 
C = 1.5 cm where the tissue response is strong, a = 0.1, 
and not delayed, r = 0. Exactly this situation matches 
the conditions under which the self-localization is most 
pronounced. When the tissue response is delayed the 
effect of the tumor existence is weakened as it is demon- 
strated in Fig. 0, where we set r = 1 min. 
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FIG. 4: The necrosis growth in living tissue containing a small 
tumor. The tissue exhibits the strong response without delay. 
The dotted lines show the position of the necrosis interface at 
different moments of time. 



IV. CLOSING REMARKS 

Concluding the paper we would like to state once more 
the obtained result. Namely, we have analyzed the necro- 
sis growth due to thermal coagulation induced by laser 
light absorption and limited by heat diffusion into the 
surrounding live tissue. The tissue is assumed to contain 
a certain tumor where the blood perfusion rate does not 
change during the action. The latter certainly concerns 
only the undamaged part of the tumor. By contrast, the 
normal tissue responds strongly to increase in the tissue 
temperature and, as a result, the blood perfusion rate 
can grow by tenfold. 

We studied in detail the necrosis formation under con- 
ditions typical for a real course of thermal therapy treat- 
ment, the duration of the action has been taken about 
5 minutes when the necrosis domain of size about or 
above 1 cm is formed. In particular, we have shown that 
if the tumor size is sufficiently large, it exceeds 1 cm, and 
the tissue response is not too delayed, the delay time does 
not exceed 1 min, then there are conditions under which 
the relative volume of the damaged normal tissue is small 



in comparison with the tumor volume after the tumor is 
coagulated totally. 
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FIG. 5: The necrosis growth in living tissue containing a large 
tumor. The tissue response to temperature increase is strong 
and not delayed. The dotted lines show the position of the 
necrosis interface at different moments of time. 
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